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I; 

ft' The generalized traveling salesman problem (GTSP) is an extension of the well-known 

traveling salesman problem. In GTSP, we are given a partition of cities into groups and we 

■ , are required to find a minimum length tour that includes exactly one city from each group. 

The aim of this paper is to present a problem reduction algorithm that deletes redundant 

vertices and edges, preserving the optimal solution. The algorithm's running time is 0(N 3 ) 

i ii in the worst case, but it is significantly faster in practice. The algorithm has reduced the 

problem size by 15-20% on average in our experiments and this has decreased the solution 

l\ , time by 10-60% for each of the considered solvers. 
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1 Introduction 

*n : 

r — I The generalized traveling salesman problem (GTSP) is defined as follows. We are given a weighted 

complete undirected graph G on N vertices and a partition V = V^UV^U. . .L)Vm of its vertices; the 
I* \ subsets Vi are called clusters. The objective is to find a minimum weight cycle containing exactly 

one vertex from each cluster. There are many publications on GTSP (see, e.g., the surveys [HE] 
and the references therein). The problem has many applications, see, e.g., 0HI]- It is NP-hard, 
since the traveling salesman problem (TSP) is its special case (when \Vi\ = 1 for each i). The 
weight of an edge xy of G is denoted dist(a;, y) and will be often called the distance between x and 

><: 

Various approaches to GTSP have been studied. There are exact algorithms such as branch- 
and-bound and branch-and-cut described in [3]. Another approach uses the fact that GTSP can 
be converted to an equivalent TSP with the same number of vertices [13l [14j [15] and then can 
be solved with some efficient TSP solver such as Concorde pQ. Heuristic GTSP algorithms have 
also been invistigated, see, e.g., [7] 151 ITU ! US ] ITS ] I3P ] l?T] . 

Different preprocessing procedures are often used for hard problems to reduce the computation 
time. There are examples of such approaches in integer and linear programming (e.g., [9l I17j) as 
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well as for the Vehicle Routing Problem [12]. In some cases preprocessing plays the key role in 
an algorithm (e.g., [5]). We introduce preprocessing procedure for GTSP. A feature of GTSP is 
that not every vertex of a problem should be visited and, thus, GTSP may contain vertices that a 
priori are not included in the optimal solution and may be removed. We have a similar situation 
with edges. 

The experimental results show that almost each GTSP instance tested in the literature can be 
reduced by the presented procedure at a very low cost and that this reduction is almost always 
beneficial for the GTSP solvers. 



2 Vertex Reduction 

Since GTSP solution covers only M vertices, up to TV — M vertices may be reduced without a 
change of the optimal solution. We present an approach to detect some of the redundant vertices 
in a reasonable time. 

Definition 1. Let C be a cluster, \C\ > 1. We say that a vertex r 6 C is redundant if for each 
pair x, y of vertices from distinct clusters different from C, there exists a vertex s S C\ {r} such 
that dist(x, s) + dist(s, y) < dist(a;, r) + dist(r, y). 

In other words, if for each path xry there exists another path xsy, s G C \ {r}, with the 
same or smaller weight, vertex r can be removed. Testing this condition for every vertex will take 
approximately 0(N 3 ■ \V\), where \V\ = N/M is the average cluster size. In the symmetric case 
of the problem there is an efficient heuristic that usually allows to reduce the preprocessing time 
significantly. 

Let us take two distinct vertices r and s in some cluster C. We can calculate the differences 
between the distances to r and s from each vertex x ^ C (A£ s = dist(x,r) — dist(x, s)) and save 
this information to a Differences Table such as Table [TJ Notice that in Table Q] we assume that 
clusters 1 and 2 have three vertices each and cluster 3 has two vertices, r belongs to the first cluster 
and it is the first vertex in the cluster, i.e., vertex s can be only the second and the third vertices 
of cluster 1. 
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Observe that a vertex r is redundant if there is no pair of vertices from different clusters such 
that the sum of differences A (see above) for these vertices is negative for every s, i.e., r is redundant 
if for every x and y there exists s £ C \ {r} such that A£ s + A^' s > 0, where x and y belong to 



distinct clusters. That is due to 



a^ + a;-' 



= dist(a;, r) — dist(ai, s) + dist(y, r) — dist(y, s) 
= dist(ai, r) + dist(r, y) — (dist(a;, s) + dist(s, y)) 

Therefore we need to check every pair of columns (col i, col 2) (except the pairs of columns 
corresponding to the same clusters) in the Differences Table T row ^ co i. If Ti tCO i 1 + Ti ]CO / 2 < 0, we 
check the second row (T2 iCO ii +T2,coi 2 )- If the result is still negative, we check the third row, etc. If 
all the rows are checked and each time we obtain a negative sum, the vertex r cannot be removed 
and the rest of the procedure may be skipped. 

Example 1. In the example above (Table [J) it is necessary to perform up to 6 tests provided in 
Table EJ 

Table 2: Vertices pairs for the example. 
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-3 




-1 


cl.2 v.2— cl.3 v.2 






4 







cl.2 v.3— cl.3 v.l 
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The only test that does not allow us to declare the vertex r redundant is in the row 3 of the 
Tabic [2] (cl.2 v.2 — cl.3 v.l) as both sums (for s = v.2 and for s — v.3) are negative. (Certainly, 
there is no need to calculate the sum for s = v.3 in rows 2, 4, and 6 in the example above, and the 
calculations may be stopped after the row 3.) 

Removing redundant vertex may cause a previously irredundant vertex to become redundant. 
Thus, it is useful to check redundancy of vertices in cyclic order until we see that, in the last cycle, 
no vertices are found to be redundant. However, in the worst case, that would lead to Q(N 2 ) 
redundancy tests. (Recall that N is the total number of vertices in GTSP.) Our computational 
experience has shown that almost all redundant vertices will be found even if we restrict ourselves 
to testing each vertex of GTSP at most twice. Thus, we assume in the rest of the paper that each 
vertex is tested at most twice for redundancy. 

2.1 Acceleration Heuristic 

In some cases it is possible to determine faster that a vertex r is not redundant. If 

min max A!j s + min max A!j s < 
x£z sec xez sec 

for some cluster Z, then r cannot be reduced. This condition means that there exist two columns 
in the Differences Table corresponding to distinct clusters and the sum of these columns maxima 
is negative. This ensures that the sum for every row of these columns is also negative. 



We can use an equivalent condition: 

min max A' ,s + min max A!}" < 
ieU K , Vj sec xev z sec 

This condition can be tested during the Difference Table generation. For each column we calculate 

the maximum value: 

vertexmax(x) = max A!) s 
sec 

Also for each cluster Z, we have 

clustermin(Z) = mmvertexmax(x) 

xez 

We define totalmin(i) — min^i clustermin(Vj); if totalmin(i) + clustermin(Vi) < for some i, 
we can conclude that vertex r is not redundant. 

In the example above, the heuristic performs just one check for V2 and V3. We have totalmin(3) = 
clustermin(V2) = —1 and clusterminiV^) — 1 and —1 + 1 > so the acceleration heuristic does 
not reduce our computations in this case. 

Another way to make the redundancy test faster is to order the rows of the Differences Table 
such that the row with the minimal number of negative values would be the first one. Notice that, 
if this row contains no negative values, it is obvious that r is redundant. 

2.2 Algorithm Complexity 

Let .Kmin and -ftT max be the minimum and the maximum number of tests (of vertices) for redundancy. 
Observe that -ftT m i n = N, since we will perform only N tests if no vertex is detected to be redundant. 
Since we have assumed that no vertex is tested more than twice for redundancy, K max = 2N — 1 . 

Now consider how many operations are required for each redundancy test (with a fixed vertex 
r). The test requires table generation and table processing. Due to the acceleration heuristic, 
table generation can be aborted already after processing of two clusters. Thus, in the best case it 
takes -E m in = (|C| — 1)(|A| + \Y\) operations where r E C, and X and Y are some other clusters. 
The average size of a cluster can be estimated as N/M (recall that M is the number of clusters). 
Therefore, in the best case each redundancy test requires approximately 

— >«(£->)('■£)-'■(£' 

operations, where N' is the current number of vertices in the problem. 

In the worst case both the table generation and the further table inspection will be completed 
normally. Table generation will take (\C\ — 1)(A' — \C\) operations. Table inspection takes about 

N' — \C\ — \V\ ) (\C\ — 1) /2 operations in the worst case, where \V\ is the average cluster size. 
Thus, we have the following number of operations per test in the worst case: 

(n'-\c\-W\Y(\c\-i) 

E max (N') « (|C| - 1) (N> |C|) + A -1 

1 ' 2 2M ' 



The total number of operations in the worst case is 

N 3 N 4 
K max ■ E max (N) « 2 ■ N ■ — = — . 



N\ 2 _ N 3 



The total operation number in the best case is 

K min -E min (N)^N-2-^j =2- u . r 

Since usually M = Q(N), the algorithm complexity changes from O(N) to 0(N 3 ). The exper- 
imental algorithm complexity is Q(N 2A ) (see Section POT) . 



3 Edge Reduction 

Definition 2. Let u, d be a pair of vertices from distinct clusters U and C respectively. Then 
the edge uv is redundant if for each vertex x G V \U \ C there exists v' € C \ {v} such that 
dist(w, u') + dist(w', x) < dist(u, v) + dist(u, x). 

Testing this condition for every edge will work for both symmetric and asymmetric cases and 
will take approximately 0(iV 3 - 1 V^|), where \V\ is the average cluster size. We introduce an algorithm 
for edge reduction for the symmetric case of the problem; it proceeds as follows. Given a vertex 
v € C, where \C\ > 1, we detect redundant edges incident with v using the following procedure: 

1. Select an arbitrary vertex v" 6 C \ {v}. 

2. Set P x = Al' v " for each vertex x <eV\C (recall that A£ s = dist(a;, r) - dist(ar, sj). 

3. Sort array P in non-decreasing order. 

4. For each cluster U j^ C and for each vertex ueJ/do the following: 

(a) S = A'r" 

(b) For each item A^ v " of the array P such that A£ ,v + 5 < check the following: if x £ U 
and A^ + A^' v < for every v' G C \ {w, w"}, the edge uw is not redundant, continue 
with the next u. 

(c) Edge uv is redundant, set dist(u,w) = oo. 

To prove that the above edge reduction algorithm works correctly, let fix some edge uv, u G U, 
v G C, U ^ C. The algorithm declares this edge redundant if the following condition holds for 
each x €" C (see |4b) : 



A v x ' v +A^ >0 or 

A v / + A v u ' v ' > for some v' G C \ {v, v"} 
This condition is equivalent to 

A v x ' v ' + A v u - V ' > for some v' £ C \ {v} 



So the algorithm declares the edge uv redundant if for each x £ V\C\U there exists v' <E C\{v} 
such that A v x ' v ' + A v u ' v ' > 0. 

dist(x, v) — dist(x, v') + dist(u, v) — dist(u, v') > and 
dist(w, v) + dist(u, x) > dist(u, v') + dist(w', x). 



Let us evaluate the algorithm's complexity. The edge reduction algorithm performs the follow- 
ing steps for every cluster C, \C\ > 1 for each v E C: 

• Array P generation. This takes 0(iV) operation. 

• Array P sorting. This takes 0(-/Vlog 2 iV) operations. 

• Edges uv testing. Each test takes 0(1) to O(7V-|0|) operations and Q(N) tests are performed. 

Thus the complexity of the entire algorithm is Q(N 2 log 2 N) in the best case, and Q(N 3 ■ \C\) 
in the worst case. 

As usually \C\ = 0(A), we may say that this algorithm's complexity varies from 0(A 2 log 2 N) 
to Q(N 3 ). The experimental algorithm complexity is 0(A 26 ) (see Section l4~Tj) . 

After the search for redundant edges has been completed, the edge reduction algorithm finds 
redundant vertices using the following observation: if after the edge reduction procedure some 
vertex has finite distance edges to at most one cluster, then this vertex can be declared redundant. 

This reduction takes 0(N 2 ) operations. 



4 Experiments 

We tested the reduction algorithms on the standard GTSP instances (see, e.g., [21 [HI HH ED] ) which 
were generated from some TSPLIB [16] instances by applying the clustering procedure of Fischetti, 
Salazar and Toth 0. The algorithms were implemented in C++ and tested on a computer with 
AMD Atlon 64 X2 Core Dual processor (3 GHz frequency). 

We have tested three reduction algorithms: the Vertex Reduction Algorithm (see Section [2]), 
the Edge Reduction Algorithm (see Section [3]), and the Combined Algorithm witch first applies 
the Vertex Reduction Algorithm and then the Edge Reduction Algorithm. 

4.1 Experimental Results 

Each test was repeated ten times. The columns of the table are as follows: 

• Instance is the instance name. The prefix number is the number of clusters of the instance; 
the suffix number is the number of vertices (before any preprocessing) . 

• R v is the number of vertices detected as redundant. 



R e is the number of edges detected as redundant. For the Combined Algorithm R e shows 
the number of redundant edges in the already reduced by the Vertex Reduction Algorithm 
problem. 

T is the preprocessing time in seconds. 

Table 3: Test results of the Reduction Algorithms. 
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13 
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41.4 


33.6 
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9dantzig42 
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0.0 


36.2 
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14.3 


24.9 


0.0 


10att48 


15 


31.3 


0.0 


41.5 
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0.0 


31.3 


25.3 


0.0 


10gr48 


18 


37.5 


0.0 


27.0 
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0.0 


37.5 


25.5 


0.0 


10hk48 
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12.5 
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0.0 


12.5 


32.3 


0.0 


llberlin52 


15 


28.8 
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25.9 
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14 


10.3 
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10.3 


26.3 
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28grl37 


10 


7.3 
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19.9 
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7.3 


17.0 
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29prl44 


19 


13.2 
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33.2 
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0.0 


13.2 


31.1 
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30chl50 


22 


14.7 


0.0 


19.9 
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14.7 


18.1 
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13.3 


0.0 


22.5 
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Vertex reduction Edge reduction Combined reduction 



Instance R v R v , % T R e , % R v T R v , %, R ei % T 

30kroB150 14 9.3 0.0 23.8 2 0.0 9.3 23.4 0.0 

31prl52 34 22.4 0.0 37.5 7 0.0 22.4 26.6 0.0 

32ul59 33 20.8 0.0 23.5 3 0.0 20.8 15.1 0.0 

35sil75 45 25.7 0.0 27.4 5 0.0 25.7 17.5 0.0 

36brgl80 97 53.9 0.0 57.9 51 0.0 53.9 16.9 0.0 

39ratl95 12 6.2 0.0 22.2 1 0.0 6.2 20.4 0.0 

40dl98 7 3.5 0.0 23.1 4 0.0 3.5 24.2 0.0 

40kroA200 16 8.0 0.0 20.3 2 0.0 8.0 20.6 0.0 

40kroB200 7 3.5 0.0 19.1 1 0.0 3.5 18.5 0.0 

41gr202 4 2.0 0.0 18.8 1 0.0 2.0 18.5 0.0 

45ts225 40 17.8 0.0 20.0 2 0.0 17.8 11.2 0.0 

45tsp225 12 5.3 0.0 20.5 2 0.0 5.3 17.1 0.0 

46pr226 12 5.3 0.0 29.6 1 0.0 5.3 28.4 0.0 

46gr229 1 0.4 0.0 22.0 0.0 0.4 21.6 0.0 

53gil262 16 6.1 0.0 21.8 3 0.0 6.1 18.9 0.0 

53pr264 11 4.2 0.0 21.5 1 0.0 4.2 20.7 0.0 

56a280 20 7.1 0.0 19.4 1 0.0 7.1 16.1 0.0 

60pr299 15 5.0 0.0 16.2 0.0 5.0 14.7 0.0 

641in318 13 4.1 0.0 20.5 2 0.0 4.1 20.8 0.0 

641inhp318 13 4.1 0.0 20.5 2 0.0 4.1 20.8 0.0 

80rd400 11 2.8 0.0 14.8 1 0.1 2.8 13.0 0.0 

84A417 43 10.3 0.0 28.3 5 0.1 10.3 22.7 0.1 

87gr431 0.0 0.0 17.2 0.3 0.0 17.2 0.3 

88pr439 10 2.3 0.0 14.7 1 0.2 2.3 15.0 0.1 

89pcb442 24 5.4 0.0 11.9 0.1 5.4 9.7 0.1 

99d493 4 0.8 0.0 17.8 1 0.2 0.8 19.4 0.2 

107att532 21 3.9 0.0 20.5 2 0.3 3.9 18.1 0.3 

107ali535 29 5.4 0.1 16.6 2 0.5 5.4 14.3 0.5 

107si535 96 17.9 0.0 26.5 9 0.3 17.9 17.9 0.1 

113pa561 147 26.2 0.1 31.3 5 0.3 26.2 22.6 0.1 

115u574 11 1.9 0.0 14.4 1 0.2 1.9 14.0 0.2 

115rat575 18 3.1 0.0 11.2 2 0.2 3.1 10.9 0.1 

131p654 88 13.5 0.1 32.6 2 0.8 13.5 28.2 0.5 

132d657 8 1.2 0.0 10.8 0.3 1.2 9.6 0.3 

134gr666 0.0 0.0 11.6 1.0 0.0 11.6 1.0 

145u724 34 4.7 0.1 10.1 3 0.5 4.7 8.8 0.4 

157rat783 25 3.2 0.0 9.8 2 0.4 3.2 8.4 0.3 

200dsjl000 8 0.8 0.1 9.6 1 2.4 0.8 9.4 1.5 

201prl002 20 2.0 0.1 9.2 2 3.0 2.0 8.7 1.6 

207sil032 85 8.2 0.2 12.1 12 1.2 8.2 10.2 0.9 

212ul060 36 3.4 0.1 14.4 1 1.7 3.4 11.2 2.0 

217vml084 241 22.2 0.6 24.0 8 2.3 22.2 8.9 1.3 

235pcbll73 11 0.9 0.1 8.2 1.5 0.9 8.2 1.3 

259dl291 48 3.7 0.2 12.4 2 2.3 3.7 9.8 1.7 
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261rll304 


19 


1.5 


0.2 


7.9 
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2.6 


1.5 


7.2 


2.0 
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23 


1.7 


0.2 


7.8 
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4.1 


1.7 


7.0 


2.9 


276nrwl379 


11 


0.8 


0.2 


7.4 
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7.1 


2.6 


280fll400 


23 


1.6 


0.9 


17.4 





6.5 


1.6 


17.5 


5.3 


287ul432 


33 


2.3 


0.2 


7.7 
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3.2 


2.3 


6.6 


2.6 


316A1577 


44 


2.8 


0.4 


10.3 


2 


5.0 


2.8 


9.2 


4.5 


331dl655 


14 


0.8 


0.2 


6.7 
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3.7 


0.8 


6.7 


3.7 


350vml748 


285 


16.3 


2.5 


19.8 
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16.3 
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5.5 


364ul817 
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0.3 


0.1 


6.2 





4.9 


0.3 


5.8 


4.5 


378rll889 


17 


0.9 


0.7 


7.3 


3 


10.9 


0.9 


6.8 


7.2 


421d2103 


8 


0.4 


0.2 


6.7 


1 


2.9 


0.4 


6.6 


2.7 


431u2152 


10 


0.5 


0.3 


5.2 





7.8 


0.5 


5.0 


6.6 


464u2319 


24 


1.0 


0.6 


3.9 





10.3 


1.0 


3.8 


9.7 


479pr2392 


33 


1.4 


0.9 


5.9 


1 


15.4 


1.4 


5.3 


13.4 


608pcb3038 


29 


1.0 


1.4 


4.7 


1 


45.4 


1.0 


4.7 


36.2 


759A3795 


21 


0.6 


4.9 


6.4 





127.2 


0.6 


6.5 


94.5 


893fnl4461 


22 


0.5 


3.4 


3.1 





80.2 


0.5 


2.9 


46.7 


1183rl5915 


28 


0.5 


7.9 


2.4 


2 


258.1 


0.5 


2.3 


114.1 


1187rl5934 


38 


0.G 


9.4 


3.0 


2 


308.3 


0.6 


2.7 


139.6 


1480pla7397 


196 


2.6 


31.5 


4.6 


1 


2147.9 


2.6 


3.6 


1001.3 


2370rlll849 


37 


0.3 


40.7 














2702usal3509 


21 


0.2 


98.7 















The results of the experiments show that the preprocessing time for the Vertex Reduction is 
negligible for all the instances up to 212ul060, i.e., for almost all TSPLIB-based GTSP instances 
used in the literature. The average percentage of detected redundant vertices for these instances is 
14%, and it is 11% for all considered instances. The experimental algorithm complexity is about 

0(N 2A ). 

The Edge Reduction is more time-consuming than the Vertex Reduction. The running time 
is negligible for all instances up to 115rat575. Note that in most of the GTSP literature, only 
instances with N < 500 are considered. The average per cent of the detected redundant edges 
for these instances is about 27%, and it is 21% for all instances in Table [3] The experimental 
algorithm's complexity is 0(N 2 - 6 ). 

4.2 Algorithms Application Results 

Certainly, one can doubt the usefulness of our reduction algorithms since they may not necessarily 
decrease the running time of GTSP solvers. Therefor, we tested the improvement of the running 
time of the following GTSP solvers: 

1. Exact algorithm (Exact) based on a transformation of GTSP to TSP [2]; the algorithm 

from 4 was not available. The algorithm that we use converts a GTSP instance with N 



vertices to a TSP instance with 37V vertices in polynomial time, solves the obtained TSP 
using the Concorde solver [I], and then converts the obtained TSP solution to GTSP solution 
also in polynomial time. 

2. Memetic algorithm from [19] (SD). A memetic algorithm (MA) is a combination of a genetic 
algorithm with local search. 

3. MA from 7 (GKK). 

4. MA from 18\ (SG). 

5. A modified version of MA from [8], the state-of-the-art GTSP memetic solver, (GK). 



Each test was repeated ten times. The columns of the tables not described in Section 1471] are 
as follows: 



• To is the initial problem solution time. 

• B is the time benefit, i.e., (To — T pr )/To, where T pr is the preprocessed problem solution 
time; it includes preprocessing time as well. 



Table 4: Time benefit for Exact. 





T , sec 


Vertices Red. 


Edge 


Red. 


Combined Reduction 


Instance 


R v , % 


B,% 


Re, % 


B,% 


R v , % 


Re, % 


B,% 


5gr21 


0.8 


38.1 


40 


48.7 


52 


38.0 


45.0 


56 


5ulysses22 


1.7 


50.0 


60 


44.3 


48 


50.0 


39.5 


79 


5gr24 


0.2 


54.2 


74 


33.1 


53 


54.1 


10.4 


81 


6fri26 


0.9 


50.0 


67 


28.7 


18 


50.0 


20.3 


74 


6bayg29 


6.0 


41.4 


19 


0.0 


59 


41.3 


33.6 


70 


10gr48 


16.1 


37.5 


57 


27.0 


2 


37.5 


25.5 


55 


10hk48 


52.7 


12.5 


16 


34.2 


6 


12.5 


32.3 


22 


lleil51 


32.8 


17.6 


37 


32.6 


17 


17.6 


28.8 


42 


14st70 


150.4 


17.1 


43 


36.5 


17 


17.1 


24.6 


50 


Average 




35.4 


45.9 


31.7 


30.2 


35.3 


28.9 


58.8 



The experiments show that the Vertex Reduction, the Edge Reduction and the Combined 
Reduction Technique significantly reduce the running time of the Exact, SD and GKK solvers. 
However, the Edge Reduction (and because of that the Combined Reduction Techique) is not that 
successful for SG (Table [7]) and the original version of GK. That is because not every algorithm 
processes infinite edges well. 



Next we show that a solver can be adjusted to work better with preprocessed instances, 
this purpose we modified GK as follows: 



For 



• The 2-opt heuristic [8] was extended with the cluster optimization. For every iteration of 
2-opt, where edges i>i«2 and V3V4 are removed, instead of replacing them with V1V3 and ^2^4 
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Table 5: Time benefit for GKK. 





T , sec 


Vertices 


Red. 


Edge 


Red. 


Combined Reduction 


Instance 


R v , % 


B,% 


Re: % 


B,% 


R v , % 


Re, % 


B, % 


89pcb442 


60.7 


5.4 


4 


11.9 


17 


5.4 


9.7 


35 


99d493 


85.2 


0.8 


14 


17.8 


19 


0.8 


19.4 


29 


107att532 


101.2 


3.9 


9 


20.5 


20 


3.9 


18.1 


20 


107ali535 


99.3 


5.4 





16.6 


47 


5.4 


14.3 


51 


107si535 


166.1 


17.9 


12 


26.5 


14 


17.9 


17.9 


41 


113pa561 


101.8 


26.2 


15 


31.3 


21 


26.2 


22.6 


47 


115u574 


103.6 


1.9 


-3 


14.4 


12 


1.9 


14.0 


28 


115rat575 


219.3 


3.1 


38 


11.2 


36 


3.1 


10.9 


45 


131p654 


165.4 


13.4 


21 


32.6 


12 


13.4 


28.2 


38 


132d657 


189.1 


1.2 


10 


10.8 


22 


1.2 


9.6 


24 


134gr666 


224.8 


0.0 


26 


11.6 


36 


0.0 


11.6 


57 


145u724 


232.9 


4.6 


25 


10.1 


29 


4.6 


8.8 


55 


157rat783 


392.7 


3.1 


1 


9.8 


16 


3.1 


8.4 


29 


200dsjl000 


898 


0.8 


6 


9.6 


52 


0.8 


9.4 


51 


Average 




6.3 


12.7 


16.8 


25.2 


6.3 


14.5 


39.3 



we replace them with v^v^ and v' 2 V4, where v[ £ cluster(v\) and v' 2 £ cluster(v 2 ) and v[ 
and v 2 are selected to minimize the solution objective value. (Here cluster(v) is the cluster 
corresponding to the vertex v. v € cluster(v).) Thereby, while the initial 2-opt heuristic 
could decline some good 2-opt if w(viV3) = oo or w(v 2 Vi) = oo, the extended 2-opt will pass 
round the infinite edges. 

• Direct 2-opt heuristic [5] is excluded from the Local Search Procedure. 

• Every time before starting the Cluster Optimization [8] we remove all vertices that cannot 
be included in the solution, i.e., if a fragment of the solution corresponds to clusters C\, C 2 
and C3 and there is no edge from C\ to v S C 2 or there is no edge from v to C3 then v can 
be excluded for the current Cluster Optimization run. 

• Since the modified Local Search Procedure is more powerful than the previous one, we 
reduced the number of solutions in a generation and the termination condition is also 
changed (now r = 0.2G + 0.03M + 8 while previously r = 0.2G + 0.05M + 10 and I cm > 
max(1.5/ max , 0.025M + 2) instead of J cur > max(1.5/ max , 0.05Af + 5), see 0). 

The modified algorithm does not reproduce exactly the results of the initial GK heuristic; it 
gives a little bit better solution quality at the cost of slightly larger running times. However, one 
can see (Table [SJ that all the Reduction Algorithms proposed in this paper influence the modified 
GK algorithm positively. 

Different reductions have different degree of success for different solvers. The Edge Reduction 
is more efficient than the Vertex Reduction for GKK and SD; in other cases the Vertex Reduction 
is more successful. For every solver except SG the Combined Technique is preferred to separate 
reductions. 
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Table 6: Time benefit for SD. 





T , sec 


Vertices Red. 


Edge 


Red. 


Combined Reduction 


Instance 


R v , % 


B,% 


Re, % 


B,% 


R v , % 


Re, % 


B,% 


157rat783 


23.6 


3.2 


11 


9.8 


5 


3.1 


8.4 


36 


200dsjl000 


100.3 


0.8 


47 


9.6 


36 


0.8 


9.4 


42 


201prl002 


54.9 


1.9 


12 


9.2 


22 


1.9 


8.7 


43 


207sil032 


21.3 


8.2 


3 


12.1 


-1 


8.2 


10.2 


24 


212ul060 


88.8 


3.3 


8 


14.4 


35 


3.3 


11.2 


42 


217vml084 


78.1 


22.2 


49 


24.0 


-2 


22.2 


8.9 


57 


235pcbll73 


107.9 


0.9 


5 


8.2 


30 


0.9 


8.2 


32 


259dl291 


169.4 


3.7 


9 


12.4 


25 


3.7 


9.8 


26 


261rll304 


140.4 


1.5 


9 


7.9 


47 


1.4 


7.2 


66 


265rll323 


132.6 


1.8 


20 


7.8 


20 


1.7 


7.0 


32 


276nrwl379 


111.5 


0.8 


4 


7.4 


22 


0.7 


7.1 


46 


Average 




4.4 


16.1 


11.2 


21.7 


4.4 


8.7 


40.5 



Preprocessing is called to reduce the solution time. On the other hand, there is no guaranty 
that the outcome of the preprocessing will be noticeable. Thus, it is important to ensure at least 
that the preprocessing time is significantly shorter than the solution time. 

Five GTSP solvers are considered in this paper. The first solver, Exact, is an exact one and, 
thus, it is clear that its time complexity is larger than 0(7V 2 - 6 ) (see Section I4TTJ) or even the upper 
bound 0(N 3 ). The time complexities of the other four solvers were estimated experimentally, i.e., 
experiments were conducted for problems of different size obtained from TSPLIB [16 and then an 
approximation for "solution time" / "instance size" dependence was found. The experimental com- 
plexity of SD is about B(N 3 ) and it is about 9(V 3 - 5 ) for GK, SG and GKK. Table M demonstrates 
the quality of our estimate for SD (here T cst imatc (N) = 6.3319 ■ 10~ 8 ■ N 3 ). 

Having the solvers time complexities, we can conclude that the preprocessing time is signifi- 
cantly smaller than the solution time for arbitrary large instances as the experimental complexity 
of preprocessing is smaller than the complexity of even the fastest of the considered solvers. 



5 Conclusion 

The GTSP reduction techniques allow one to significantly decrease the problem complexity at a 
very low cost. Experiments show that the Combined Reduction is often the most powerful among 
the presented algorithms and takes even less time than the single Edge Reduction. While the 
Vertex Reduction yields very natural problems and is successful with every considered solver, the 
Edge Reduction changes some edge weights to infinity values and, thus, not every solver benefits 
from it. However, in this paper, it is shown that a solver can be modified to process such problems 
well. 

In this paper we consider the symmetric case only, i.e., dist(x,y) = dist(y, x) for every pair of 
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Table 7: Time benefit for SG. 





To, sec 


Vertices Red. 


Instance 


Rv, % 


B,% 


84fi417.gtsp 


4.5 


10.3 


12 


87gr431.gtsp 


8.3 


0.0 


6 


88pr439.gtsp 


10.2 


2.3 


-3 


89pcb442.gtsp 


11.5 


5.4 





99d493.gtsp 


20.0 


0.8 


7 


107att532.gtsp 


25.1 


3.9 


11 


107si535.gtsp 


16.9 


17.9 


34 


107ali535.gtsp 


29.1 


5.4 


20 


113pa561.gtsp 


14.5 


26.2 


31 


Average 




8.0 


13 



vertices x and y. Other vertex and edge reduction algorithms that can be immediately derived 
from Definitions [T] and [2] exist for the asymmetric case, and their time complexity is 0(N 3 ). Recall 
that N is the total number of problem vertices. 

Acknowledgement We would like to thank Larry Snyder and John Silbcrholz for kindly providing 
the source codes of SD and SG, respectively. 
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Tabic 9: SD work time estimation. 
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60pr299 

641in318 

641inhp318 

80rd400 

84H417 

87gr431 

88pr439 

89pcb442 

107si535 

113pa561 

115u574 

115rat575 

131p654 

134gr666 

145u724 

157rat783 

200dsjl000 

201prl002 

207sil032 

212ul060 

217vml084 

235pcbll73 

259dl291 

261rll304 

265rll323 

276nrwl379 
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